Asymptotic properties of wall-induced chaotic mixing in point vortex pairs. 
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The purpose of this work is to analyze the flow due to a potential point vortex pair in the vicinity 
of a symmetry line (or "wall"), in order to understand why the presence of the wall, even far from 
the vortices, accelerates fluid mixing around the vortex pair. An asymptotic analysis, in the limit 
of large distances to the wall, allows to approximate the wall effect as a constant translation of 
the vortex pair parallel to the wall, plus a straining flow which induces a natural blinking vortex 
mechanism with period half the rotation period. A Melnikov analysis of lagrangian particles, in the 
frame translating and rotating with the vortices, shows that a homoclinic bifurcation indeed occurs, 
so that the various separatrices located near the vortex pair (and rotating with it) do not survive 
when a wall is present. The thickness of the resulting inner stochastic layer is estimated by using the 
separatrix map method, and is shown to scale like the inverse of the squared distance to the wall. 
This estimation provides a lower-bound to the numerical thickness measured from either Poincare 
sections or simulations of lagrangian particles transported by the exact potential velocity field in the 
laboratory frame. In addition, it is shown that the outer homoclinic cycle, separating the vortices 
from the external (open) flow, is also perturbed from inside by the rotation of the vortex pair. As 
a consequence, a stochastic layer is shown to exist also in the vicinity of this cycle, allowing fluid 
exchange between the vortices and the outer flow. However, the thickness of this outer stochastic 
zone is observed to be much smaller than the one of the inner stochastic zone near vortices, as soon 
as the distance to the wall is large enough. 
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I. INTRODUCTION 

Vortical flows have remarkable transport properties 
which have been investigated in various contexts in the 
past decades. Indeed, vortices are often referred to as 
"violent" or "singular" structures which induce large 
local velocity gradients which in turn significantly in- 
fluence the transport of scalar or vector fields. For ex- 
ample, steady two-dimensional isolated fixed vortices, 
whether pointwise or with a viscous core, have been 
shown to significantly influence the transport of pas- 
sive scalars (Flohr & Vassilicos (H) or magnetic fields 
(Bajer [2]). If, in addition, the vortex oscillates, then 
the velocity field is unsteady, and this unsteadiness is 
likely to influence scalar mixing through some chaotic 
advection phenomenon (Aref Q |H, Beigie et al. Q). 
One of the simplest unsteady vortex of this kind is the 
periodically translating isolated vortex, which can be 
thought of as a variation of the blinking vortex, and 
which has been shown to induce fast mixing in chaotic 
regions (Wohnas & Vassilicos @). In the case of vor- 
tex pairs the flow is also unsteady, since both vortices 
rotate around each other (provided the total strength 
is non-zero), but this unsteadiness is not sufficient to 
induce chaotic mixing. Indeed, in the frame rotat- 
ing with the vortex pair the flow is steady, and the 
motion of fluid points is an autonomous dynamical 
system with only 1 degree-of-freedom: their motion 
cannot be chaotic. 

If, in addition, the two vortices move in the vicin- 
ity of a flat wall, then elementary vortex dynamics 
shows that both the distance between the vortices and 



their rotation velocity are no longer constant. Fluid 
point trajectories are therefore much more complex. 
It has been shown that chaotic advection indeed oc- 
curs in the case of leapfrogging vortex pairs (Pentck 
et al. which can be thought of as a vortex pair 
near a wall simulated by using the mirror method. In 
the present paper we will investigate this phenomenon 
analytically to show that chaos persists when the vor- 
tices are far from the wall, and to quantify this effect 
by calculating the thickness of the resulting stochastic 
layers. 

We will use an asymptotic method to show that the 
basic ingredient of chaos in the work by Rom-Kedar 
et al. [8J , namely the artificial unsteady straining flow 
applied to a counter-rotating vortex pair, is also at 
work in the case of a co-rotating vortex pair near a 
wall. These authors show that counter-rotating vor- 
tex pairs can induce chaotic advection, provided they 
are superposed to an oscillating straining flow. They 
observe that in the absence of the oscillating strain- 
ing flow "the pair carries a constant body of fluid", 
whereas the oscillating straining flow forces the fluid 
to be "entrained and detrained from the neighbour- 
hood of the vortices". 

The goal of the present paper is to show that, for 
co-rotating vortex pairs, a wall-induced straining flow 
can make the two vortices oscillate to-and-fro with a 
period half the natural rotation period of the pair, and 
that this variation of the blinking vortex phenomenon 
triggers chaotic mixing. This result will be derived by 
using an asymptotic approach, in the limit do -C Lq, 
where do is the order-of-magnitude of the distance be- 
tween vortices, and Lq is the distance to the wall. This 
"inner" asymptotic velocity will approximate the flow 
in the vicinity of vortices, i.e. at distance O(do) from 
them. Wall-induced chaotic mixing will then be ana- 
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lyzed by considering the breaking of the corresponding 
"inner" homoclinic cycles. 

In addition, it will be shown that vortex rotation 
also perturbs, from inside, the outer homoclinic cycle 
separating the vortex pair from the external open flow. 
This effect will be analysed by means of an "outer" 
asymptotic velocity field (in the same limit do ^ Lq), 
designed to approximate the flow at distance Lq from 
vortices. 

The inner asymptotic model is depicted in section 
HI! Chaotic advection and mixing in this flow are in- 
vestigated and discussed in section IIIII We then fo- 
cus on fluid exchange with the outer open flow by 
using the outer approximation (section IIV|) . Because 
do <C Lo, the motion of vortices will be assumed to be 
regular throughout the paper (see Acheson Q). 
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FIG. 1: Sketch of the vortex pair {A\, A2), together with 
the two mirror vortices {Ai, A' 2 ). 
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where (x\,yi) are the coordinates of point A\, and 
(^2,2/2) are the coordinates of point A^. The latter 
sum is the wall effect, and will be next simplified in 
the limit where the wall is far from the vortices (see 
for example Appendix A of Rom-Kedar et al. 
Elementary vortex dynamics shows that: 
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yi = constant = Lq , 

in all cases, with Zk = Xk+iyk- By setting yk = yi+yk, 
with y~k = O(do) <C Lq, and expanding ([5]) we get: 

x I =x a I (l + 0(e 2 )), 



with x j 



T/2ttLq, and e = cIq/Lq <C 1. This shows 
that, to leading order, / moves under the effect of 
a single mirror vortex with strength — 2T. If non- 
dimensionalized by Oo and do, the x velocity of I 
reads: 



VLodo 



= 2e + 0(e 3 ), 



(3) 



so that the vortex pair moves at constant speed up 
to order 3. In order to exploit further the fact that 
do Lo, we consider positions (x, y) located in the 
vicinity of the vortex pair, and normalize the variables 
as: x* = (x - x°j)/d , y* = (y - yi)/d , and x* = 
(xi-x°j)/do, y* = (y l -yi)/d , where x°j = Tt/2irL . 
Also, we set t* = Qot and define the non-dimensional 
streamfunction in the frame translating at the velocity 
of /: 



r(x*,y*,t*) 



ip-x j(y- yj) 



(4) 



II. INNER ASYMPTOTIC VELOCITY FIELD 
NEAR VORTICES 



Then, by expanding the streamfunction in the limit 
e — > 0, we get (removing the stars): 



We consider a pair of point vortices {Ai,Ai) with 
equal strength T moving in the plane (O, x, y). In the 
absence of any solid boundary both vortices rotate 
around their centre point / with an angular velocity 
Do = T/4irdQ, where 2<io = l^i^bl is the (constant) 
distance between the vortices in this case. In the pres- 
ence of a wall, located at y — 0, the streamfunction is 
(inviscid fluid): 



i>(x,y,t) = ^ -Log [(a 



(y - Vif 



+ -(x 2 -y 2 ) + 0(e 3 ). 



(5) 



Terms of order 0(e) have cancelled-out with 
expected. Note that this flow does not correspond to 
a uniform flow — at infinity (in contrast with Eq. 
^ ) , as it is only an approximation of the streamfunc- 
tion in the vicinity of the vortices. 
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A. Approximate vortex motion 

We therefore observe that the effect of the wall can 
be approximated by a stretching flow with axes e x ±e y 
(where e x and e y are the unit vectors attached to axes 
x and y respectively). Note that this flow is very 
close to the one investigated by Carton et al. [10( 
and Maze et al. The line A X A 2 joining the two 

vortices is stretched and compressed periodically dur- 
ing the motion. It will be shown later that the pe- 
riod of this stretching/compression is half the period 
of the vortices without wall. The dynamics of vortex 
A\ = (xi,yi) in this simplified flow is: 



±1 = -2 
Vi = 2 



yi - m 

X\ — x 2 



4r 2 



-{Vi ~ V2), 



— (x 1 -x 2 ), 



(6) 
(7) 



with 2r(t) = \ AiA 2 \. By writing the motion equation 
of A 2 = (x 2 ,y 2 ) we are led to: (i\ + x 2 )/2 = = 
(z/i +2/2)/2. This shows that, to this order, the centre 
point / is indeed fixed, in agreement with equation 
©■ The latter equation was also expected: indeed, 
as noticed in the introduction, yj is known to remain 
constant whatever the distance to the wall. The mo- 
tion of the vortices can then be investigated further 
by setting x\(t) = r(t) cos9(t), yi(t) = r(t) sin#(i), 
so that equations © and ([?]>. with x 2 — —x\ and 
y% = —yi lead to: 



-e r sin 29, 
e 2 cos 29. 
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(8) 
(9) 



Looking for an asymptotic solution in powers of e, we 
get (setting r(0) = 1 + £ 2 /2 and 9(0) = 0): 



"(t) = I + — cos2i, 



9{t) 



t - 



■ sin 2t, 



(10) 
(11) 



and the vortex coordinates in the reference frame 
translating with / then read: 



xi(t) = cost + e 2 



cost 



sini 



cos2i + sinisin2t (12) 



(t) = sint + e 2 cos 2t- cost sin 2t (13) 



yi( 



We therefore observe that the wall, when far from the 
vortex pair, creates an oscillation of the distance be- 
tween the two vortices with period 7r (half the period 
of the isolated vortex pair). Also, the rotation velocity 
of the vortices around I is affected by a perturbation 
with period tt. This periodic forcing corresponds to 
the effect of the straining flow £ 2 (— y, — x) appearing 
in the velocity field, which stretches and compresses 
the vortex pair twice during each rotation (see Fig. 




FIG. 2: Sketch of the stretching flow induced by the wall in 
the reference frame translating with /. This flow stretches 
and compresses the segment A\Ai twice during each com- 
plete rotation of the vortex pair. 



In the following lines we write x = (x,y), and 
the vortex position is also written in vector form: 
Xi(t) — x^ +e 2 x^ 2 \ with the components of the var- 
ious vectors given by (|T2")) - (rr3"|) . Then the streamfunc- 
tion in the frame translating with /, written in powers 
of e, reads: 

2 

ip a (S,t) = ^2 tpo (2 -Xi(i)) +e 2 ^ 2 {x), (14) 

i=l 

with ipo(x) = — Log(a; 2 +y 2 ) and tp 2 (x) = (x 2 - y 2 )/2. 
It can be expanded as: 



+ e 2 il) 2 (x). 



(15) 



This streamfunction corresponds to the flow induced 
by the two vortices and the wall, in the reference frame 
translating with /, and is valid only in the vicinity of 
the vortex pair A\A 2 . 



B. Flow in the rotating frame 

It is very convenient to investigate the motion of 
lagrangian particles in the reference frame translating 
with / and rotating at speed fio; that is the angu- 
lar velocity of the vortex pair in the absence of a wall. 
Let ex , eV denote an orthogonal basis attached to this 
frame, with X, Y the corresponding coordinates. As- 
suming ex — S x and eV = e y at t = we have: 
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x = X cos t — Y sin t and y — X sin t + Y cos t. Then, 
in this reference frame, and with these coordinates, 
the streamfunction reads 

MX, Y, t) = </> a (x(X, Y, t),y{X, Y, t), t) + ^(X 2 +Y 2 ), 

where the last term corresponds to the entrainment 
velocity of the rotating frame. By using (fl~5l) we ob- 
tain, after some algebra: 

MX,Y,t) =Ao(X,Y)+e 2 A2(X,Y,t), (16) 

where 

lM*,Y,t) = -Log|Z 2 - 1| 2 + i|Z| 2 , (17) 

and 

yj r2 (X,Y,t) = - 2 [(Y 2 -X 2 + l)cos2i 



+ 4AY"sin2i] + - [(X 2 -Y 2 )cos2t~2XYsm2t] , 

where Z = X + iY. The leading order streamfunction 
corresponds to the steady flow in the rotating frame 
in the absence of wall, the e 2 terms are the wall effect. 



III. WALL-INDUCED CHAOTIC 
ADVECTION NEAR VORTICES 



the occurrence of simple zeros in the Melnikov func- 
tion of the separatrix E,: 

/CO 
[Vvo (£(*)), ^V-a (&(*), * + to)} dt, (20) 
-co 

where to is the starting time of the Poincare section of 
the perturbed tracer dynamics, and [., .] is the Poisson 
bracket. By choosing <fo(0) on the "middle-point" of 
Si (i.e. the intersection between the separatrix and 
OY for Ei and E 2 , and OX, X ^ 0, for E ), then the 
coordinates of qi.(t) are either odd or even functions, 
and by using the various symmetries of the velocity 
field, one can get rid of the integral of odd functions 
of t in (HO]). We are led to Mj(i ) = Q-i sin(2t ), where 
the amplitudes of the Melnikov functions oti are purely 
numerical constants. They are calculated by solving 
the (fi(i)'s numerically on the three separatrices (in 
the upper half-plane for £1,2 and in the right half- 
plane for So): a n ~ —0.58, ot\ ~ —0.89, a-i ~ 7.31. 
We therefore conclude that the two homoclinic cycles 
and the two homoclinic trajectories of Fig. [3] do not 
survive when a wall is present. This is confirmed by 
the Poincare sections shown in Fig. SJ A stochastic 
zone indeed exists in the vicinity of the separatrices, 
and grows with e. It will be called "inner" stochastic 
zone in the following, as it exists in the very vicinity 
of the vortices (i.e. at distance O(do) from them). 
The typical scale of this zone, and its boundaries, are 
investigated in the next section. 



A. Splitting of the inner homoclinic cycle 

In order to investigate the occurrence of chaotic ad- 
vection in the flow (|16[) we consider fluid points, the 
dynamics of which reads: 

dX dipr _ di> r _> 

~dT = W ex - dx eY > (19) 

where X(t) — (X, Y) is the position of the fluid point. 
This is a perturbed hamiltonian system, the phase 
portrait of which is, to leading order, the same as the 
one of Fig. 2a of Angilella [12}, which is reproduced 
here (Fig. [3]) for clarity. This widely used phase por- 
trait has two homoclinic cycles (Si U S2 and its sym- 
metric in the lower half-plane), and two homoclinic 
trajectories (So and its symmetric in the left half- 
plane). Let qi(t) be a solution of the leading-order 
fluid points dynamics: 

dqi _ dlp r o _ dlp r o _ 
~tt-^Y ex ~^X eY > 

with initial condition on a separatrix E,; (i = 0, 1, 2). 
Then the occurrence of chaos in the fluid point dy- 
namics under the effect of the e 2 terms is related to 
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FIG. 3: Streamlines of the leading-order flow in the rotat- 
ing frame, in the vicinity of the vortices (inner asymptotic 
field). 
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FIG. 4: Poincare sections of 20 particles, together with the 
theoretical boundary of the inner stochastic zone (black 
thick line). Graph (a): e = 0.05, graph (b): e = 0.2. 



B. Characterization of the inner stochastic zone 

In order to quantify mixing in the vicinity of the vor- 
tices we derive a semi-analytic expression for the equa- 
tion of the border of the stochastic layer. The thick- 
ness of this zone grows with e, and can be calculated 
from the above Melnikov functions (Chirikov fl3T |. 
Weiss & Knobloch [14] . Rom-Kedar [la], Kuznetsov 
& Zaslavsky [l6j], Trueba & Baltanas [l7| . Balasuriya 
[18]). We build the separatrix map of £, for i g {1, 2} 



(Fig. [S]) by considering a fluid point X(t) running un- 
der the effect of the perturbed flow ([TUl) . This point 
crosses the axis IY at discrete times t n . In addition, 
following Weiss & Knobloch [TJ], we define the set 
S composed of the hyperbolic points A and B, and 
of segments forming two symmetric stars near A and 
B (Fig. [5|). These stars have two branches parallel 
to IX, and two other branches lying between Si and 
£2. The detailed shape of the S is of no importance to 
construct the separatrix map, provided the trajectory 
of X(t) remains close to the separatrices. In addition 
we call r n the time at which the trajectory of X(t) in- 
tersects S just prior time t n , that is: r n < t n < r n +i. 

Let H n denote the value of the unperturbed stream- 
function ?/vo at X(r n ). The variation of "energy" H n 
between two consecutive crossings of S is: 



AH n = H, 



n+l' 



ip r o(X(t)),Tp r2 (X(t),t) 



If the trajectory is close enough to a separatrix 
(Chirikov [13]) , say E^, then we have X(t) ~ %{t— t n ), 
where <fj, parametrizing Ej, has been introduced in 
the previous section, and is such that <fi±(0) be- 
longs to the vertical axis IY. Also, by noticing that 
T n <C t„ <C T n +i, we obtain: 



(21) 



where M{ are the Melnikov functions of separatri- 
ces Ej calculated in the above section. The time lag 
t n+ i — t n can be approximated by the half-period of 
the unperturbed orbit [l3| ip r o = i?„+i, so that we set: 



tn+l — tn + T(H n+ i) /2. 



(22) 



Equations (|2"Tj) and ([2^]) approximate the separatrix 
map in the vicinity of the homoclinic cycle E1UE2. To 
our knowledge, no exact analytical expression for the 
period T(H) is known. However, following Kuznetsov 
& Zaslavsky pj], one can notice that the dynamics of 
fluid points is very slow in the vicinity of point A, so 
that an approximate expression of T(H) can be found 
by expanding the streamfunction in the vicinity of A. 
We then obtain hyperbolic streamlines for this sim- 
plified dynamics, and T(H)/2 is approximated as the 
time required for a fluid point running on an hyper- 
bolic streamline to pass over A. We are led to (see 



Appendix : 



T(H) ~ -PLag\H A - H\ 



(23) 



with /3 = 4/V5. 

The amplitude of the rate of change of the time 
interval d{t n+ i — t n )/dt n can then be calculated from 
equations and (|2"Tj) : 



Max 



dt 



n+l 



dt n 



s 2 l3\ ai \/\H A -H n+1 \ 



(24) 



dt. 
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FIG. 5: Sketch of the separatrix map used to calculate the 
boundary of the inner stochastic layer. 



By writing that this quantity is equal to unit y on the 
border of the stochastic layer (Rom-Kedar [la ]) we 
obtain the energy H s corresponding to this boundary: 



\H.-H A \ =P\a l \e 2 



(25) 



The boundary of the stochastic layer around £2 can 
then be obtained by setting i = 2 in the above for- 
mula. The curve i/vo = H s has been plotted on the 
numerical Poincare sections of Fig. @] for comparison 
(thick black line), and we observe that, in spite of the 
rough approximations of the model, this curve corre- 
sponds to the border of the stochastic zone with an 
acceptable accuracy. 

Result (|25|) can also be used to determine the dis- 
tance re between the border of the stochastic layer 
and point C, which can be thought of as the half- 
thickness of the stochastic layer at C . Indeed, by writ- 
ing that -ip r o(Xc + rc<=V) = H s . and using a first-order 
Taylor expansion (assuming that tq is small): 



r c 



(3a2 2 
7 



(26) 



where 7 = dtl> r o/dY(C) ~ 1.67. The thickness of 
the stochastic layer of £2 at C is therefore: Ac = 
2rc — 15. 8e 2 . This formula is compared to measure- 
ments of the thickness Ac obtained from numerical 
computations (Fig. [6]). Three kinds of computations 
are used here: Poincare sections, particle clouds com- 
puted by using the asymptotic velocity in the rotat- 
ing frame, and particle clouds computed by using the 
exact potential velocity |T]) in the laboratory frame. 
Point C has been chosen for these measurements be- 
cause the thickness is well-defined there. We observe 



that the order-of-magnitude of Ac is well predicted 
by the theory, especially if e is small. For larger e dis- 
crepancies of about 100% appear, and these errors are 
known to be due to resonances occurring in the fluid 
point dynamics (Luo & Han Trueba & Baltanas 
|r7j |) and to regular islands embedded in the stochastic 
layer. In spite of these discrepancies we can conclude 
that the theoretical prediction obtained from the sep- 
aratrix map gives an acceptable lower bound to the 
thickness of the stochastic zone. 




FIG. 6: Thickness of the inner stochastic zone at point C 
of E2, versus e = do/Lo, obtained from the separatrix map 
(equation (|26l) . thick solid line) and estimated from numer- 
ical computations: Poincare sections (empty circles), par- 
ticle cloud using the asymptotic velocity in the rotating 
frame (squares), particle cloud using the exact potential 
velocity in the laboratory frame (filled circles). The dashed 
line is the theoretical thickness of the outer stochastic zone 
near the separatrix S1S2 discussed in section HVl 

The typical size of the stochastic layer can also be 
estimated in the vicinity of the hyperbolic points A 
(and B) , by writing that iproi^A + Ta^x) = H s , 
where ta is the distance from A to the external bor- 
der of the stochastic zone, towards ex. To obtain ta 
we then use a second-order Taylor expansion (assum- 
ing that ta is small). Indeed, because A is a stag- 
nation point of the unperturbed flow, the gradient of 
?/vo vanished there, and we are left with: 



r A 



/3a2 
5 



1/2 



(27) 



The typical size of the stochastic layer is therefore 
larger in the vicinity of A (and B) than anywhere else, 
in the limit where e is small. This is a rather general 
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result which is a consequence of the fact that A and 
B are stagnation points [lj| . Here also we expect (J57J 
to provide only a lower bound of the actual , since 
resonances might occur and make the thickness even 
bigger. 



C. Comparison with transport by the exact 
potential velocity field 

The asymptotic model presented in the above sec- 
tions shows that some chaotic mixing must occur in 
the vicinity of a vortex pair as soon as a wall is present 
in the surroundings of the two vortices. In order to 
check the reliability of these calculations, which are 
valid only if the wall is far from the vortices, we have 
performed some calculations of fluid points transport 
by using the exact potential velocity field ([1]) (set di- 
mensionless by using Oo and do), in the laboratory 
frame. 

Fig. [7] shows two particle clouds corresponding to 
Lq = 10do and 5do (that is e — 0.1 and 0.2, respec- 
tively). Because our goal is to visualize fluid point 
mixing in the stochastic zone, the initial cloud is taken 
to be a small spot of size 0.1 x 0.1 centered at I 
(which is initially located at (0,0)). The initial po- 
sitions of the vortices are (x±,y±) = (1 + e 2 /2,l/e) 
and (x2,U2) = (—1 — £ 2 /2,l/e). The final time of 
the three simulations is t = 50ir, and the final loca- 
tion of the centre point is very close to the asymptotic 
value (xi,yi) « (2et, Lq) in both cases. We observe 
that some particles collect on a structure which is very 
close to the separatrices Sj described above, as a con- 
sequence of the existence of a stochastic zone in the 
vicinity of these separatrices. The thickness of the 
stochastic zone is of the order of the theoretical one, 
as shown in Fig. [5] (filled circles). 



IV. MIXING WITH THE OUTER OPEN 
FLOW 

The results of the above sections show that chaotic 
advection can occur in the vicinity of the vortices, due 
to the presence of the wall which induces a periodic 
unsteadiness in the vortex motion. However, the inner 
asymptotic model (valid only at distance O(do) from 
vortices) does not bring any information about fluid 
mixing with the "external flow" which carries fluid ele- 
ments from the far field. Fig.JSJa) shows the complete 
flow (solid lines), obtained from the exact potential 
solution (JTJ), in the frame translating with vortices. 
Both the real and mirror vortices are shown. Lengths 
and velocities have been normalized by Lq and T /Lq 
respectively, which are relevant scales for the exter- 
nal flow, and which will be used to derive an outer 
asymptotic velocity field below. The vortex system is 
separated from the rest of the flow by a homoclinic 
cycle which can be thought of as the boundary be- 
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FIG. 7: Particle clouds at t — 50n, computed in the lab- 
oratory frame by using the exact potential velocity field, 
for L = Wd (e = 0.1) (a), and L = 5d (e = 0.2). Solid 
lines are streamlines. 



tween closed and open streamlines. The breaking of 
this homoclinic cycle is obviously an important mech- 
anism for fluid mixing between the vortex system and 
the fluid outside (Pentek el al. [7]). This homoclinic 
cycle will be called the "outer cycle" in the following, 
as opposed to the so-called "inner cycle" £2 U Si. 

Because the flow is two-dimensional, the only can- 
didate for the splitting of this outer cycle is the ro- 
tation of the vortices around each other, which in- 
duces an unsteady perturbation to the homoclinic cy- 
cle Si 5*2- This bifurcation is investigated in the next 
paragraphs. 



<s 



r, r, — r, — r, in this frame, is: 





FIG. 8: Graph (a): streamlines in the frame translating 
with the vortices, when e = 0.3. Both the real (y > 0) and 
mirror (y < 0) vortices are shown ; the wall corresponds 
to y = 0. Solid lines: exact four-vortex potential flow 
JT|). Dashed lines: outer asymptotic streamfunction (|29[) . 
Graph (b): leading-order flow (e = 0) in the same frame. 
Si and S2 are saddle points located at (±\/3, 0). The 
thick line is the homoclinic cycle S\Si. (Lengths have 
been rescaled by Lq, stars have been removed.) 



A. Asymptotic analysis of the outer homoclinic 
cycle 

As noticed above, Lq and T being the relevant 
scales of the external flow, we non-dimensionalize 
the streamfunction ([T]) by using Lq for lengths and 
T/AttLq for velocities. Also, we consider the refer- 
ence frame translating at the velocity Xj — T/2itLq 
of a vortex with strength 2r, placed at distance Lq 
from the wall. The non-dimensional streamfunction 
of the flow induced by the four vortices with strengths 



T/Att ' 



(28) 



where ij)(x,y,t) is the streamfunction of the flow ob- 
served in the laboratory frame, and stars indicate 
non-dimensional variables. Still assuming that e = 
do/L <C 1, and using Eqs. (fT2")l and P^|) for the vor- 
tex motion, the streamfunction can be expanded as: 

2t 

ip*(X*) ~ V>o*(^*) +e 2 V'2c(^*)cos— y 



+ £ 2 V^(X.)sin^ + 0( £ 4 ), (29) 

where -0o*(^*,y*) = 21og[(x 2 +y 2 + I + 2y*)/(a; 2 +y 2 + 
I — 2y*)] — 2y* is the streamfunction of a simple dipole 
centered at (0,0) (i.e. single vortex plus its mirror 
vortex, see Fig. Efb)), and the e 2 terms manifest the 
fact that we have a vortex pair with strength r + T, 
instead of a single vortex with strength 2r, so that the 
flow is unsteady. The streamfunctions tp2c and tp2s are 
steady, rational functions of x» and only. 

With the time scale used here, the rotation of the 
vortex pair is fast, so that the flow (|2T)|) is a rapidly 
perturbed dipolar flow. Fig. OJa) shows this stream- 
function when e = 0.3 (dashed line), together with the 
exact four-vortex streamfunction already discussed 
(solid line). The two flows are very close, except in 
the vicinity of both pairs, as expected. The 0(e 3 ) 
terms in (|29p are identically zero (for symmetry rea- 
sons), and this is why the agreement is good even if 
e is not very small here. For e < 0.2 dashed lines 
and solid lines (not shown here) are indistinguishable, 
except in the very vicinity of the vortices. 

Particle motion can therefore be approximated by a 
perturbed hamiltonian system, with hamiltonian ^0* . 
The Melnikov function attached to the cycle S1S2 
reads (removing the stars): 



M ext (t ) 



[Vto (${*)), iMtft*))] cos (2 



t + t 



)dt 



[Vo(<7(t)), iMtft*))] sin(2^^)dt, (30) 



where q{t) runs on the unperturbed cycle, and is cho- 
sen such that <f(0) coincides with point C of Fig. IHJb). 
This Melnikov function approximates (to order e 2 ) 
the "energy" difference ipo(X(T2)) — ipo(X(Ti)), where 
X(t) is the trajectory of a fluid point advected by the 
perturbed flow, with a trajectory very close the sep- 
aratrix S\S2- The times t\ and t-i correspond to the 
particle passing near the saddle points S\ and 8%. Be- 
cause our perturbation is rapid, the amplitude of the 
Melnikov function will now depend on e (Gelfreich 
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[20]). The integrands in the above integrals are un- 
known analytically, they are however localized around 
zero and decay exponentially. We have chosen to ap- 
proximate them with a mixed polynomial/gaussian 
expression of the form P(t)/ cosh 2 (5 1), where P is a 
fitted polynomial and 8 is a fitted constant. The in- 
tegrals can then be obtained analytically, and we are 
led to M ext (t ) = A(e)sm(2t /e 2 ), with: 



A(e, 



(31) 



where 6 and 7 are purely numerical constants depend- 
ing on the fitting procedure: 7 18, 5 « 5. Note that 
this formula depends on the form chosen to approxi- 
mate our integrands. This is therefore a numerical re- 
sult, rather than a truely analytical one, being under- 
stood that other fittings might not significantly affect 
the numerical values of the amplitude of the Melnikov 
function. Because A(e) > for all e, and M ext (to) al- 
ways has simple zeros when t$ varies, portions of fluid 
will be exchanged between the interior of the cycle 
S1S2, and the fluid outside. However, the amplitude is 
so small as soon as e < 0.2 (i.e. A(0.2) < 7 10~ 5 ), that 
mixing through the cycle Si 62 should be hardly visible 
for such small e's. In contrast, when e increases above 
0.2, A(e) quickly reaches values of order unity. By 
constructing the separatrix map of Si S2 , one can de- 
rive the theoretical thickness A of the stochastic layer 
surrounding the outer homoclinic cycle. In the vicin- 
ity of point C of Fig.HKb), we are led to A c « 0.5A(s). 
This quantity has been plotted in Fig. [5] (dashed line). 
We indeed observe that the thickness is exponentially 
small when e < 0.2 (Lq > 5do), and quickly reaches 
values of order unity as soon as e is larger than about 
0.25. [Note that, even though the diagram suggests 
a sharp increase of Ac, there is no critical e here: 
the Melnikov function has simple zeros for all e] Fig. 
IH1 suggests that the vortex pair might be apparently 
separated from the outer open flow when Lq > 5do, 
since the whole vortex system will be surrounded by 
an extremely thin stochastic zone. Also, the thickness 
of this stochastic zone should quickly reach a value of 
order unity as soon as Lq < 3do- These two points are 
checked numerically in the next paragraph. 



B. Numerical verification 




(d) 
e = 0.3125 




Computations of particle clouds in the laboratory 
frame, using the exact four-vortex potential solution 
((T|) , have been performed. Particles have been injected 
at t = near the vortices (in box A of Fig. HJb)). 
Also, to check whether portions of fluid could enter 
within the cycle S1S2, we have injected another set 
of particles outside the cycle S1S2 (in box B of Fig. 
EJb)). These simulations confirm the significant differ- 
ence between the thickness of the inner cycle and the 
one of the outer cycle, as suggested by the theoretical 
curves of Fig. [SJ 



FIG. 9: The role of the inner and outer separatrices in 
chaotic mixing. The four graphs show particle clouds 
in the laboratory frame at t = 1007r, obtained by using 
the exact four-vortex potential solution. Particles marked 
with black dots have been injected near the upper vortices 
at t = (in box A of Fig.[8jb)). Those marked with empty 
circles have been injected outside the cycle S1S2, ahead 
of vortices (in box B of Fig. MP))- Graph (a): e ~ 0.2 
(Lo/do = 5). Graph (b): e ~ 0.25 (L /d = 4). Graph 
(c): e ~ 0.286 (Lo/do = 3.5). Graph (d): e ~ 0.3125 
(Lo/do — 3.2). Solid lines are streamlines. 
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Indeed, Fig. M, a ) shows the two particle clouds at 
t = 1007T, when Lq/cIq = 5 (e — 0.2). Units are the 
same as in the simulations of section Hill that is: d 
for lengths and Slo for times. Particles marked with 
black dots are those which have been injected in box A 
(they are also indicated with a (A)), and those marked 
with empty circles have been injected in box B (they 
are also indicated with a (B)). Even though these 
computations have been conducted in the laboratory 
frame, the homoclinic cycles appearing in the rotating 
frame (section Mil) and in the translating frame (sec- 
tion |IVX|, naturaly appear here. Particles injected 
near vortices wander in the vicinity of the inner cy- 
cle T,2 U Ei investigated in section Hill whereas those 
injected ahead of the vortices move in the vicinity of 
the outer cycle S11S2. Some particles penetrate within 
this cycle, in agreement with the fact that the Mel- 
nikov function has simple zeros. When Lo/do = 4 
(Fig. Ejfb)) the thickness of the inner stochastic zone 
around the cycle £2 is larger, whereas the thickness 
of the outer stochastic zone around the cycle S1S2 
does not appear to be very different from case (a), in 
agreement with the theoretical thickness calculated in 
the previous section (dashed line of Fig. When 
Lo/do = 3.5 (Fig. [9jc)) the thickness of the outer 
stochastic zone is slightly larger than in case (b), but 
the two clouds are still separated. For Lo/do = 3.2 
(e = 0.3125, Fig. (Hid)) the two particle clouds mix sig- 
nificantly, in agreement with the fact that the ampli- 
tude of the Melnikov function is of order unity there. 
However, the thickness of the outer stochastic zone 
cannot be measured from these graphs and compared 
to the theoretical value of Ac- this is due to the fact 
that particles injected in box B now mix with those 
injected in box A. These simulations confirm, how- 
ever, that the outer stochastic zone remains negligible 
compared to the inner one as long as Lo/do is larger 
than about 4, so that the vortex pair might exchange 
very little fluid with the outer flow in this case. 



V. CONCLUSION 

We have analyzed how point vortex pairs, which 
are closed flow structures in the absence of walls, ex- 
change portions of fluid with their surroundings as 
soon as a wall, even far from the vortices, is present. 
This is not a new result, since Pentek et al. 0] had 
shown that chaotic advection is present in leapfrog- 
ging vortex pairs, but the asymptotic approach used 
here shows that the mechanism at work is a combina- 
tion of the stretching flow induced by the wall and of 
the unsteadiness due to vortex rotation. Our work has 
therefore some common features with the one of Rom- 
Kedar et al. Q, even though the flow is different. Rel- 
evant quantities, like the thickness of stochastic layers, 
have been calculated by using Melnikov functions, as 
done in previous works for other flows (e.g . Kuznetsov 
& Zaslavsky [l6j|, Trueba & Baltanas |T7|). 



The breaking of the inner separatrices, located in 
the very vicinity of the vortex pair and rotating with 
them, induces a stochastic layer the thickness of which 
scales like the squared inverse distance to the wall. 
The rotation of vortices around each other can also 
affect the outer homoclinic cycle S1S2 separating the 
vortex system and the rest of the flow domain. The 
thickness of this outer stochastic layer has been shown 
to be negligible as long as e < 0.25, so that the vortex 
pairs might be apparently separated from the open 
flow zone when Lq > 4do- However, this thickness 
rapidly jumps to unity for larger e's. As a conse- 
quence, significant portions of fluid coming from in- 
finity can penetrate within this cycle and mix with 
the fluid located near vortices. 

One could argue that the values of e considered in 
the various computations presented here are not very 
small, so that the comparison with the asymptotic re- 
sults might be done with care. However, the asymp- 
totic analyses show that the leading-order wall effect 
scales like e 2 . The approximation is especially accu- 
rate in the case of the outer separatrix (section IIVI) 
since the velocity truncation error scales like e 4 there. 

The generalization to vortices with unequal 
strengths can be done by using the same asymptotic 
approach when the wall is far from the vortices, pro- 
vided the sum of the two strengths is non-zero. Here 
also we obtain that the center of vorticity / (strength- 
weighted averaged position) moves at constant speed 
up to order 3 and that the wall induces a stretch- 
ing flow, in the frame translating with I, which trig- 
gers a blinking- vortex mechanism. One can check that 
the various separatrices do not survive when a wall is 
present in this case also. 

Finally, it should be noted that the "wall" analyzed 
in this work is in fact a symmetry line, rather than a 
wall, as the inviscid fluid considered here slips on solid 
boundaries. The generalization of these results to 
real fluids, using the Navier-Stokes equations together 
with no-slip boundary conditions, must be done with 
care. Indeed, it is well known that vortices, even with 
a large Reynolds number, induce counter-rotating sec- 
ondary vortices when a wall is present. This vortic- 
ity, which is initially located inside the wall boundary 
layer, can contaminate the flow and interact with the 
primary vortices, leading to a complex vortex system, 
very different from the ideal one investigated here. 
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Appendix A: Calculation of the period near the 
separatrix 

In the vicinity of saddle point A the unperturbed 
streamfunction can be expanded as: 

H = i>ro(X) ~H A + ~(J? - X A ).V 2 iP r0 .(X - X A ) 



where a is an arbitrary initial height of the fluid point, 
which has very little effect on T(H) when the trajec- 
tory is close to the separatrix. The vertical velocity 
is Vy = -(5/2)(X - X A ), where X is given by ([AT]) , 
and we get: 



H A - 



Y 2 



-\x-x A y 



(Al) 



It corresponds to open hyperbolic trajectories like the 
one shown in Fig. HUT a). For such a trajectory the 
time required to pass near point A is very large if 
SH = \H — H A \ is small, and this time is a good ap- 
proximation of the half-period of trajectories around 
the homoclinic cycle [16J. We therefore write 



T(H)/2 - 



~ a dY 
Vy 



TiH) = -^Lo, 



2SH 1/2 



a+ {ASH + G 2 ) 1 / 2 ' 



In the limit SH — > 0, keeping a fixed, we obtain Eq. 
(ESD- Fig.HUTb) shows a plot of Eq. gjft together with 
a numerical measurement of T(H), and we observe 
that the agreement is excellent. 
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FIG. 10: Graph (a): linearized flow in the vicinity of stag- 
nation point A. Graph (b): comparison between the theo- 
retical period of trajectories near the separatrix (Eq. 
dashed line) and the numerical one (dots). 



